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ABSTRACT 

We obtain the projected mass distributions for two Subaru/Suprime-Cam fields in the southwest 
region (r ^ 60.') of the Coma cluster (z = 0.0236) by weak lensing analysis and detect eight subclump 
candidates. We quantify the contribution of background large-scale structure (LSS) on the projected 
mass distributions using SDSS multi-bands and photometric data, under the assumption of mass-to- 
light ratio for field galaxies. We find that one of eight subclump candidates, which is not associated 
with any member galaxies, is significantly affected by LSS lensing. The mean projected mass for seven 

subclumps extracted from the main cluster potential is (M^™"^) = (5.06 ± 1.30) x W^H'^Mq after 
a LSS correction. A tangential distortion profile over an ensemble of subclumps is well described by a 
truncated singular-isothermal sphere model and a truncated NEW model. A typical truncated radius 
of subclumps, rt ~ 35 /i~^kpc, is derived without assuming any relations between mass and light for 
member galaxies. The radius coincides well with the tidal radius, ~ 42 ft.~^kpc, of the gravitational 
force of the main cluster. Taking into account the incompleteness of data area, a projection effect and 
spurious lensing peaks, it is expected that mass of cluster substructures account for 19 percent of the 
virial mass, with 13 percent statistical error. The mass fraction of cluster substructures is in rough 
agreement with numerical simulations. 

Subject headings: galaxies: clusters: individual: Coma Cluster (A1656) - gravitational lensing: weak 
- X-rays: galaxies: clusters 



1. INTRODUCTION 

The cold dark matter (CDM) paradigm predicts the 
presence of numerous substructures in dark halos on 
any scale, because less massive objects form earlier and 
become more massive through mergers. Indeed, high- 
resolution N-body simulations have shown an assemble 
history that subhalos continually fall into larger halos. 
When interior subhalos penetrate into a central region 
of a massive parent halo, subhalos are disrupted by its 
strong tidal field. As a result, the original subhalo mass 
is reduced by the tidal effect and becomes a part of 
smoothed component of the parent halo (e.g De Lucia 
et al. 2004; Gao et al. 2004). Therefore, a study of 
subhalo properties is of vital importance to understand 
an assemble history in halo environments. Furthermore, 
a statistical study for subhalos, such as mass function, 
would offer a powerful test of the CDM model on scales 
less than several Mpc. 

Gravitational lensing analysis on background galaxies 
is the unique technique to map out mass distributions of 
any object, such as galaxies and clusters, regardless of 
the dynamical state. It therefore enables us to explore 
substructures in primary halos and to measure directly 
their masses. Indeed, galaxy-galaxy lensing studies in 
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clusters (e.g. Natarajan & Springel 2004; Natarajan, 
De Lucia & Springel 2007) revealed cluster substructures 
and measured their mass functions under the assumption 
of a scaling relation between mass and light. However, 
a technique requiring no assumption of mass-to-light re- 
lation is of paramount importance, because subhalo size 
in a strong tidal field depends on its orbit parameters 
as well. Furthermore, we cannot rule out a possibility 
that gaseous galaxies in a gaseous environment of galaxy 
clusters are offset from subhalo centers because of ram 
pressure. Even a slight offset prevents us from measur- 
ing subhalo masses accurately, because a mis-centering 
of tangential distortion profile causes large error in mass 
estimations especially within inner regions (e.g Yang et al 
2006; Johnston et al. 2007). It is therefore of importance 
to explore subhalos and measure their masses based on a 
weak mass reconstruction technique independent of any 
mass-to-light scaling relations. 

As demonstrated by Okabe & Umetsu (2008), a sys- 
tematic weak lensing study on seven mering clusters in 
the range of z ~ 0.055 — 0.28 is capable of discovering 
massive substructures associated with cluster major ma- 
jors. However, the limit of angular resolution of recon- 
structed mass distributions, within the redshift range, 
makes it difficult to discover less massive substructures 
associated with cluster galaxies. On the other hand, it 
would be easier to detect less massive substructures in 
lower redshift clusters in spite of the weakness of their 
lensing signal, because the number of available source 
galaxies increases thanks to larger apparent size of ob- 
jects at lower redshift. Therefore, weak lensing study of 
low-redshift clusters will provides us with a good oppor- 
tunity to detect and measure smaller subhalo masses in 
clusters. As the first step, we select Coma cluster for the 
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target to measure subhalo masses by weak lensing anal- 
ysis alone. The redshift of Coma cluster is 0.0236 and is 
known as one of the most massive clusters near us. We 
analyze archival Subaru/Suprime-Cam data (Miyazaki 
et al. 2002) to measure subhalo masses found in pro- 
jected mass distributions as well as cluster virial mass, 
and calculate the mass fraction of substructures. We also 
investigate lensing from background large-scale structure 
(LSS) in a quantitative way, using SDSS multi-bands and 
photometric data. 

The outline of this paper is as follows. We briefly de- 
scribe the data analysis in f}2] and measure the three- 
dimensional mass enclosed within a spherical region of 
a given radius using a tangential shear profile in Sj3] SI 
represents projected distributions of mass and member 
galaxies, quantifies false lensing peaks and estimate back- 
ground lensing effects on the weak lensing mass recon- 
struction. In ^ we measure the two-dimensional masses 
for subclumps with and without LSS lensing correction. 
In ^j6l we fit a tangential shear profile over an ensemble 
of subclump candidates and obtain the typical truncated 
radius and mass of subhalos. f}7] is devoted to the dis- 
cussion. Throughout the paper we adopt cosmology pa- 
rameters flmo — 0.27 and ft\ = 0.73. At the redshift of 
Coma cluster 1'. = 20.0/i^^kpc. 

2. DATA ANALYSIS 

We retrieved two i?c image data (Yoshida et al. 2008) 
from the Subaru archival data (SMOKA^). Pointings of 
imaging data are the central region oi r ^ 30'. from cD 
galaxy NGC4874 and the outskirts region of r 30- 60f. 
They cover the southwest part of this cluster. The 
data were reduced by the same imaging process of using 
standard p ipeline reduction software for Supri me-Cam, 
SDFRED (lYagi et al. 20021 : lOuchi et al. 2001 . as de- 
scribed in Okabe & Umetsu (2008). Astrometry and the 
photometric calibration were conducted using the Sloan 
Digital Sky Survey (SDSS) data catalog. The exposure 
times are 42 and 16 minutes for the central and outskirt 
regions, respectively. 

Our weak lensing analysis was done using the IM- 
CAT package provided by N. Kaiser(Kaiser, Squires & 
Broadhurst 1995^). We use the same pipeline as Ok- 
abe et al. (2009) with some modifications followed by 
Erben et al. (2001) (also see Okabe & Umetsu 2008). 
In the pipeline, we first measure the image ellipticity, 
Ca, from the weighted quadrupole moments of the sur- 
face brightness of each object and then correct the PSF 
anisotropy as e'„ = - PSi {Psm*)^^e'^*, where Pap is 
the smear polarizablity tensor and the asterisk denotes 
the stellar objects. We fit the stellar anisotropy kernel 
iPsm*)afi^^* with the second-order bi-polynomials func- 
tion in several subimages whose sizes are determined 
based on the typical coherent scale of the measured 
PSF anisotropy pattern. We finally obtain the reduced 
shear ga = 7a/(l - k) = {Pg)~le'i^ using the pre-seeing 
shear polarizablity tensor Pg . We adopt the scalar value 
{Pg)ap = Tr[Fg](5Q/3/2, following the technique described 
in Erben (2001). 

^ http:/ /smoka. nao.ac.jp/index.jsp 

® http: / /www.ifa.hawaii/kaiser /IMCAT 



We ran the pipeline for each imaging data and ob- 
tained the shear catalogue of source galaxies whose mag- 
nitude ranges are 20 — 25 ABmag and half-light-radius 
are -I- cr(r^) < r/j < 10 pixel, where and cr(»'^) 
are the mean and la error for stellar objects, respec- 
tively. Here the upper limit of magnitude is determined 
by the outskirt data of short exposure time, although 
faint galaxies in the range of 25 — 26AB mag are usable 
in the data of central region. Since apparent sizes of 
unlensed galaxies, mainly cluster members, are large in 
general, our source galaxy selection efficiently excludes 
member galaxies which dilutes lensing strengths. The 
number density of source galaxies is ~ 23 arcmin"^. 

3. CLUSTER MASS MEASUREMENT 

We measure a tangential shear component, = 
—gi cos 2ip — g2 sin 2ip and the 45 degree rotated com- 
ponent, gx = —gi sin2(^ -|- 52 cos2(^, with respect to the 
cluster center, where ip is the position angle in counter 
clockwise direction from the first coordinate. Then, 
the profiles of shear components = (g+^g^) are ob- 
tained from the weighted azimuthal average of the dis- 
tortion components of source galaxies as {gi3)(9n) = 
'}2i'^g,i90,il ^i'^g,i with a statistical weight Ug^i — 
l/(crg j -l-a^), where subscripts 'n' and 'i' denote the n-th 
radial bin and i-th source object. We adopt the soften- 
ing constant a = 0.4 which is a typical value of the mean 
rms (Tg over source galaxies. 

There are two cD galaxies (NGC 4874 and NGC 4889) 
in the central region of Coma cluster. We adopt the cen- 
ter of Coma cluster as NGC 4874 because a number of 
luminous galaxies are concentrated around NGC 4874 in 
our optical image, and the peak of X-ray surface bright- 
ness is close to NGC 4874 (see also Figure [TU)) . The shear 
profile covers the range of 4f — 60.' with 5 bins. It cor- 
responds to the first bin of Kubo et al. (2007). We fit 
the shear profile with the universal profile proposed by 
Navarro, Frenk & White (1996; hereafter NEW profile) 
and a singular isothermal sphere (SIS) halo model. We 
assume that the redshift of source background galaxies 
is (zg) = 1. An uncertainty of source redshift in mass 
estimates is negligible because the lens distance ratio, 
Dis/Ds, at such a low redshift cluster weakly depends 
on source redshifts. 

The NEW mass model is described by two parameters 
of the three-dimensional mass Mnfw(< ''a) and the halo 
concentration ca — r\/rs, where Vg is a scale radius and 
TA is a radius at which the mean density is A times the 
critical mass density, Pcr(-z), at the cluster redshift. The 
density profile of NEW mass model is expressed in the 
form of 

PNFW(?') = -j—, rj^^ -I y2- (1) 

[r/rs)(l + r/rsY 

The three-dimensional cluster mass for NEW mass model 
is obtained by 

A/nfw(< ''a) = 5 to(ca), (2) 

Ca 

with 

m(a;) =log(l -t- x) - — - — . (3) 
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Fig. 1. — Profiles of tangential shear component (top panel), 
and the 45 degree rotated component (bottom panel), gy. 
The solid and dashed lines are the best-fit NFW and SIS models, 
respectively 



For reference with other works, resuhs of mass and con- 
centration parameter within radii for A = 2500, 500 & 
200 and virial overdensity A — Avir — 98 (Nakamura & 
Suto 1997), are hsted in Table [TJ The density profile for 
SIS mass model is given by 



Psisir) 



27rGr2 



(4) 



where the one-dimensional velocity dispersion, a^, is a 
parameter. 

The resulting parameters are summarized in Table [TJ 
Since the covering area of the data is small ^60'., the 
NFW mass is not constrained well. This is why the mass 
Ma determined by fitting the tangential shear is sensitive 
to the tangential distortion at (Okabe et al. 2009). 
Since the best-fit virial radius is ~ 98.'5, we require data 
of wider region to measure the cluster mass accurately. 
We note that the is quite small because the number of 
background galaxies is scarce and the intrinsic ellipticity 
noise is large. The covering area of our data is only 
^ 16% within 60'.. If the data covers whole area, the error 
would improve ^ 2.5 times and the would become 
close to 1. The NFW virial mass changes —0.5% and 
-1-0.8% if the mean source redshift is changed to (zg) = 
1.2 and 0.8, respectively. 

We also perform a fitting taking into account a large- 
scale structure (LSS) lensing effect. The estimation of 
LSS effect on lensing signal will be described in detail in 
^4.41 The best-fit parameters are summarized in Table [21 
The LSS effect is not significant on cluster mass estimate. 

4. PROJECTED DISTRIBUTIONS OF MASS AND 
GALAXIES 



4.1. Mass, Number Density and Luminosity Maps 

We reconstruct the lensing convergence field, k, from 
the shear field, using the Kaiser & Squires inversion 
method (Kaiser & Squires 1993), following Okabe & 
Umetsu (2008). In the map making, we pixelize the shear 
data into a pixel grid using a Gaussian smoothing ker- 
nel Wg{9) oc exp[—9'^/9g] and a statistical weight Ug^i 
(^. The shear field at a given position (6) is obtained 

where the weak limit ga ~ 7q is assumed. We em- 
ploy the smoothing FWHM = V4 In 26lg = 2. '00. The 
error variance for the smoothed shear is given as o'|(^) — 

(E. y^giO - 0^?ulyg,^)/ (E. ^^(^ " 0^)Ug,^f ■ Wc havC 

used (ffa,^ff/J,J) = il/2)aliS^^Sfj with ,5^^ and Sfj being 
the Kronecker's delta. In the linear map-making process, 
the pixelized shear field is weighted by the inverse of the 
variance. 

The resulting E-mode and B-mode map of lensing fields 
are shown in panels B and C in Figures [2| and [3| which 
cover the central region (26.' x 26.') and the outskirts 
(21.' X 21.'), respectively. Contours are spaced in a unit of 
1(7 reconstructed errors. As seen in panels A and B, we 
find eight candidates of mass clumps whose significance 
is over 3cr level. The panel A in Figures [2] and [3| show the 
Subaru i?c-band images overlaid with contours of recon- 
structed mass distributions. We labeled subclumps as 
shown in panel A. The significance levels of mass clumps 
are lower than those of other clusters at redshift range 
z - 0.055 - 0.28 (Okabe & Umetsu 2008). The B-mode 
map (panel C) in the central region shows 2 clumps over 
3cr close to clump candidates 4 and 2 (3.5cr and 3.9a). 

We retrieve the SDSS DR7 catalogue (Abazajian et 
al. 2009) from SDSS CasJobs site'^ in order to inves- 
tigate member galaxy distributions. We select bright 
member galaxies by criteria of i' < 19 ABmag and 
\{g' - i') - (-0.05i' + 2.04) | < 0.14, where we use psf- 
Mag for magnitude and modelcolor for color. We convert 
from apparent to absolute magnitudes by using the k- 
correction for early-type galaxies under the assumption 
that all member galaxies are located at a single cluster 
redshift. The galaxy luminosity and density projected 
distributions are obtained using the same kernel of weak 
lensing mass reconstruction. The overall mass distribu- 
tion appears to be similar to the galaxy luminosity and 
density distributions. In particular, the 6 out of 8 mass 
candidates, but for the clumps 3 and 5, host bright galax- 
ies. At the clump candidate 3, groups of faint member 
galaxies were known (Conselice & Gallagher 1999), while 
any galaxy group are not found at the candidate 5 region. 
We list the luminous galaxy associated with each candi- 
date in Table[3l We do not always detect mass structures 
for all known groups or luminous galaxies. There are 
three possibilities for this. First, the large-scale-structure 
(LSS) lensing effect prevents us from detecting lensing 
signals. Second, dark matter halos associated with al- 
most member galaxies are less massive than the detection 
limit (3cr), iyXirOpl^SK) ~ ixW^h^^MQ. Third, dark 
halos lost their mass by the tidal force of the main cluster 
and then are smoothly distributed within the smoothing 
scale of mass reconstruction. 

http://casjobs.sdss.org/ 
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TABLE 1 

Best-fit Mass Models without a LSS lensing correction 
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Note. — Col. (1): Over-densities A = Avir, 200, 500, & 2500. Col. (2): Best-fit three-dimensional cluster mass for NFW mass model 
within TA at which the mean interior density is A times the critical mass density, pci{z), at the cluster redshift (eq. |2)l. Col. (3): Best-fit 
concentration parameter ca. Col. (4): Radius ta corresponding to the mass Mnfw(< ''a)- Col. (5): Angular size of radius r/\ Col. (6): 
Best-fit velocity dispersion for SIS model (eq. |4]l. The x'^/d.o.f is the chi-square for best-fits and the degree-of-freedom. 



TABLE 2 

Best-Fit Mass Models with a LSS lensing correction (see [[3]& i|4.4li 
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Note. — See the caption of Table [T] 



4.2. Bootstrap Re-sampling Mass Reconstructions 

We run 1000 bootstrap simulations for making mass 
reconstruction in order to investigate the realization of 
mass clumps. In each reconstruction, we generate a boot- 
strap data-set by choosing randomly galaxies, with re- 
placement, from the original shear catalogue and then 
identify mass clumps whose significance level is more 
than 3(7. Figure |4] shows the resulting distributions of 
histogram of the appearance of mass peaks These dis- 
tributions are well associated with mass clumps. The 
radii at which 68% of the centroid positions contain are 
0'.8 — 3.'. Therefore, the detected lensing peaks are real- 
ized well in the shear catalogue. 

4.3. Monte- Carlo Realizations 

We next construct a noise map, Krms, from 1000 Monte- 
Carlo realizations, following Miyazaki et al. (2007). The 
position and shear components of background galaxy cat- 
alogue are randomly shuffled in each realization. A mass 
map for a new background catalogue is reconstructed 
by applying the same procedure as making the original 
K maps. We estimate the rms noise in each pixel and 
make the noise maps, Krms for the central region and the 



outskirts. Noise maps are not changed even if we ran- 
domly choose half of catalogue. The significance maps, 
v — k/ki-ois, are obtained by dividing the original k maps 
by the Kims map. The resulting v and Urms maps for both 
E- and B-modes are shown in Figure [5] The variation 
along the left and top edges of Urms map in the central 
region is smaller than the other region. This is why fewer 
galaxies exist around the boundary of the optical image. 
The significance maps, v, are consistent with the origi- 
nal E- and B-mode maps. The significance for subclump 
candidates, i', is also consistent with the original S/N 
ratio (Table [3]). 

Arc positions of E-mode and B-mode peaks correlated 
? In the central region, two B-mode peaks whose sig- 
nificance level is above 3a are appeared close to E-mode 
peaks. We calculate the probability, Peb, as function of 
the distance between E- and B-mode peaks appeared in 
Monte-Carlo realizations. The following result does not 
change even when we use half of realization data. Since 
the appearance probability is proportional to the area 
size, we also compute the probability, Prnd, that E- and 
B-mode peaks are randomly and independently located 
in each pixel. Figure [6] shows the appearance probabili- 
ties, Peb, for the central region and the outskirts, which 
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roughly agrees with the probability of white noise case. 
A slight excess of the ratio Peb / -Prnd is found in the range 
of 6* < lOf, but the probability is quite small. Since they 
are not significant, we cannot identify fake E-mode peaks 
using the distance from B-mode peaks. The probabilities 
of large distance > 20.' is smaller than the unity, because 
few peaks are appeared around the edge. Indeed, the 
appearance probability in one pixel within 2.' width from 
the boundary is about one-thirds of that in the rest re- 
gion. The probability of spurious lensing peaks will be 
eval uated considering the large-scale structure lensing ef- 
fect EH 

4.4. Projection Effects 

Since Coma cluster is quite close to us, we cannot 
rule out a possibility that lensing signals by background 
structures significantly contribute to the observed ones. 
In this subsection, we quantify the projection effect by 
background structures on local convergence peaks ap- 
peared in the mass maps, based on the observational 
data, rather than a theory. In this paper, we use the 
shear catalogue derived from one pass-band data alone, 
which makes quite difficult to measure the contributions 
of background structures on lensing signals. 

The SDSS DR7 data (Abazajian et al. 2009), on the 
other hand, allows us to quantify the contribution, be- 
cause a huge multi-band data with photometric redshifts 
are available. We retrieved the data in the region of 
10° X 10° (190° < RA < 200° and 23° < DEC < 33°). 
Since there is no candidate for galaxy clusters or groups 
at higher redshift in the Subaru data field by visual 
checks, at least, we expect to ignore contributions of 
background clusters/groups in our data field. We quan- 
tify the projection effect by field galaxies with photo- 
metric catalogue under the assumption of mass-to-light 
scaling relations (Guzik & Seljak 2002). First, we se- 
lect galaxy catalogue by r' < 21, and Zph — zi > Sz = 
o'v.maxi^ + zi)/c ~ 0.01, taking into account the uncer- 
tainty of the photometric redshift due to line-of-sight 
velocities of member galaxies. Here, Zph is the photo- 
metric redshift of each galaxy, c is the light velocity and 
o'v,max — 3000 is the maximum of the linc-of-sight ve- 
locity (Rines et al. 2003). The following results do not 
change even when we choose the redshift range of 0.5Sz 
and 26 z, because a relative contribution of low redshift 
galaxies in the lensing signal is quite small. The result- 
ing galaxy catalogue has a peak at Zph ^ 0.5 — 0.6 in 
the histogram of photometric redshifts for faint galax- 
ies (20 < r' < 21) in DR6 data (Oyaizu et al. 2008), 
If a spectroscopic data of a galaxy is available, we uti- 
lize a spectroscopic redshift instead of photometric one. 
Next, we calculate the multi-band luminosities (u'g'r'i'z') 
within a radius of 1.93 deg from each position of galaxy 
in faint shear catalogue, corresponding to 30 Mpc at 
z — 0.5, in order to consider contributions from unknown 
clusters/groups around z = 0.5, because the two-halo 
term in the tangential shear measurements (Seljak 2000; 
Mandelbaum et al. 2005) are dominated around a few 
tens of Mpc (Johnston et al. 2007). Third, we calculate 
individual galaxy masses from the multi-band luminosi- 
ties (u'g'r'i'z') assuming the mass-to-light ratios obtained 
by SDSS bands (Guzik & Seljak 2002). The assumed 
mass-to-light ratio is in agreement with results in other 
bands (Hoekstra et al. 2005). Since we adopt the mass 



and luminosity scaling relation for a galaxy, a mass of 
an overdensity region at which a distribution of galaxies 
is concentrated might be overestimated. Finally, we as- 
sume the mass-concentration relation (Duffy et al. 2008) 
to estimate NEW shear signals at each galaxy position in 
background shear catalogue, and add them all up. The 
luminosity scaling relations in multi-bands are comple- 
mentary with each other to calibrate the lensing signals 
from background large-scale structures. If the assumed 
mass-luminosity relation is adequate, the reduced shears, 
ga, estimated in each band should coincide with those in 
other bands. The resulting shears in u'g'i' bands are sys- 
tematically inconsistent with all other bands, while the 
shears in r'z' bands have a tight correlation. We here- 
after adopt g^ct^^'^ = {ga ^ + ga as a, model of lensing 
signals from background large-scale structures. 

We reconstruct the lensing convergence fields with the 
same kernel of mass maps ( M.ip using LSS contributed 
shears alone. Here, we do not add intrinsic shape noises 
for galaxies as well as shears from main cluster, in order 
to investigate the LSS lensing effect only. Figure [7] shows 
the resulting E- and B-mode maps. The signal-to- noise 
ratio for background LSS convergence field is at ~ 2ct 
level. In subclump candidates 5 and 8, peaks of ^ la are 
found in the LSS field, while no galaxy concentration in 
SDSS catalogue is found there. It indicates a possibility 
that an appearance of two clumps in the mass maps is 
biased by a LSS lensing effect. 

4.5. Probability of Spurious Lensing Peaks 

We next investigate a probability to detect spurious 
lensing peaks by a composition of the LSS and main 
cluster lensing signals and intrinsic shapes. We con- 
sider shears composed of ga = ffi™'*'"'' + Sa"^^' + ea^*'\ 
where gi™^'"' is a best-fit NFW model in ^ and Cq"*'' 
is an intrinsic shape. We produce the intrinsic elliptic- 
ities with a Gaussian distribution whose meari value is 
|ea 1 = and the standard deviation is |<5eQ | = 0.28, 
corresponding to observed shear distributions. We then 
reconstruct mass maps and identify the mass peaks above 
3a level, as the same procedures ( §4.ip . We repeat this 

1000 times. The probability, Pspu^\ to detect spurious 
lensing peak within a smoothing scale of each mass clump 
candidate except main cluster center is summarized in 
Tabic [3l The probabilities for clump candidates 5 and 

8 are P^^^f} = 12.2% and P^^^^l = 1.8%, respectively. 
They are much higher than those for the other clump 
candidates Pspuf^ < 0.4%. We perform the same steps 
without LSS lensing signals {ga — ffi™^'"'' + ei'"'''). The 
probabihties for candidates 5 and 8 are pfpur^s^ = 0.7% 

and Pspul^s^ = 0.9%, respectively. The probability to 
count spurious peaks in candidate 5 region becomes sig- 
nificantly higher by LSS lensing effect. We also compute 
the averaged probability in the region excluding clump 
candidates in the Monte Carlo simulation taking into ac- 
count the main cluster and intrinsic noises. The resulting 

averaged probabilities are (Pspuf"') = 0.30 and 0.22% in 
central and outer regions, respectively. Here, we assume 
that the realization for spurious peaks is Possionian. The 
probabilities for candidate 5 and 8 regions are still higher 
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RA [deg] 

Fig. 2. — The panel A is Subaru i?c-band image of the central 26f X 26'. cluster region. Two cD galaxies (NGC 4874 and NGC 4889) 
are located around the northeast boundary. Overlaid arc contours of the reconstructed projected mass distribution, spaced in a unit of la 
reconstruction error (Sk = 0.018). The Gaussian FWHM is 2.'00. The identified subclumps are labeled in the panel A. The panels B and 
C are the lensing k (E-mode) and B-mode fields. The panels D and E are cluster luminosity and density distributions in SDSS i'-band 
smoothed to the same angular resolution of the mass map, respectively. Seven clump candidates are found in the central region. 




Fig. 3. — The panel A is Subaru Rc-band image of the outskirts 21'. x 21.' cluster region. Overlaid are contours of the reconstructed 
projected mass distribution, spaced in a unit of la reconstruction error (<5k = 0.020). The Gaussian FWHM is 2f00. The identified 
subclumps are labeled in the panel A. The panels B and C arc the lensing k (E-mode) and B-mode fields. The clump 8 is far ~ 58^ from 
cD galaxy NGC4874. The panels D and E are cluster luminosity and density distributions in SDSS i'-band smoothed to the same angular 
resolution of the mass map, respectively. One clump candidate is found in the outskirts. 
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Fig. 4. — The histogram maps of centroid of mass peaks, whose significance level is above 3cr, appeared in 1000 bootstrap re-sampling 
mass-reconstructions. Contours arc the same as Figures [2] and |3] Left and right panels are the central and outskirts regions, respectively. 
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0.06 



0.04 



0.02 




E— /B— mode peak distance : [arcmin] 

Fig. 6. — Top: The appearance probabilities as a function of 
distance between E- and B-mode peaks for the central region and 
the outskirts. Peb is the appearance probability obtained by 1000 
Monte-Carlo realizations. Pmd is the probability that E- and B- 
mode peaks are randomly and independently appeared. Middle: 
The ratio of PEB/Pmd for the central region. Bottom: The ratio 
of ^'Es/^rnd for the outskirts. 

even except LSS effects. It might be due to distributions 
for sheared galaxies. We define the bias for a prefer- 
ence to detect a clump in weak mass reconstruction as 

^map = ^^s''pur^V(^spur^°^)- As Summarized in Tabled the 
bias for candidates 5 and 8 are significantly higher than 
those for others. The appearances for candidates 5 and 
8 in the reconstructed mass maps are likely to be due to 
background LSS lensing distortions, at least partially. In 
the next two sections, we will quantify this more accu- 
rately based on two complementary mass measurements 
using shears, because each pixel in the convergence field 
is correlated by smoothing kernel. 

5. PROJECTED MASS MEASUREMENTS FOR SUBCLUMPS 

We measure projected mass of subclump candidates 
(labeled in Figures [2] and [3]) which are identified above 
the 3ct significant level in the mass map. The projected 
masses, M^^(< 9), are estimated by the so-called Cc- 
statistics (Clowe et al. 2000) which is modified version 
of original one (Fahlman et al. 1994), 



(5) 



where. 



d\n9 (7+)(^) 



1 



-k{< 9) -k 



< 



d\ne (7+)((?) 
? < ^b2). (6) 



Here, is a given radius and, 0bi and 0b2 are the inner 
and outer radii of subtracted background region. The 
critical surface density Ecr = c^/(47rG)(Ds£'is/-Di) is 
given by the angular diameter distances to cluster (-Di), 
to source {Ds) and between cluster and source {D\s). The 
Mq^ is a mo del- independent mass estimation. 

We first obtain central positions of mass clumps 
by peak-finding algorithms and then redistribute them 
within 0.'2 over 500 Monte Carlo simulations to take into 
account uncertainties of central positions. We choice the 
central position at which the signal-to-nose ratio of Cc 
measurement is at the maximum (Table [3]). The back- 
ground region is in the annulus of 40 — lOOkpc surround- 
ing candidates so that we extract the substructure mass 
embedded in the cluster main potential. If the clus- 
ter potential is uniform within lOOkpc, it is a good 
mass estimate of subclumps. The following results are 
not changed by choosing the background region. As de- 
scribed in eq. ([H]), Cc mass measurement is computed 
by integrating the measured tangential shears outside 
a given radius 9. Since the available number for back- 
ground galaxies in the outer annulus is larger than those 
in the inner one, it enables us to plot Mq^ profile for each 
mass clump. 

We also calculate LSS-corrected projected masses, 
in terms of (7+) {9) — (7^'"^^^) (0), where 7^^^^ is 
the LSS shears obtained in i)4.4l We do not consider the 
intrinsic shape noise. We note that M^'^"'^'^-' is a model- 
dependent mass because we assumed the mass-to-light 

ratio in a calculation of LSS shears, 7!'"^^^ Figured] 
shows the LSS-uncorrected and corrected projected mass 
profiles, respectively. The values of Af^^ are saturated on 
the outer radius, which indicates that the mass density 
of clumps is quite low on the outer radius. We esti- 
mate two-dimensional masses for each clump, M2D and 
M^"^ , from the saturated values with a covariance ma- 
trix, because each bin is correlated with each other. We 
list the resulting masses in Table [3l The LSS- uncor- 
rected and corrected masses are in good agreement with 
each other, but for candidate 5. The LSS-corrected mass 
for candidate 5 is about half of that before a correction. 
Therefore, the background LSS effect significantly con- 
tributes on lensing signals for candidate 5. 

The mean values are (M^™"^) = (4-84 ± 1.41) x 
IO12/1-IM0 for all candidates and {m!^°"^) = (5.06 ± 
1.30) X W^hr'^MQ but for candidate 5. They are at 
the order of the mass of cD galaxy halos, m'^^\9) — 
3.0 X 10i2^-iMQ(cr.„/400km/s)2(r/40kpc), where cr„ is 
the velocity dispersion and we employ the velocity dis- 
persion of cD galaxy ^ 200 — 400km/s (e.g. Smith et al. 
2000). Although we looked into the SDSS photometric 
and spectroscopic data, there is no correlation among the 
projected mass, and velocity dispersion and luminosity 
of galaxies which are located in each subclump. 

6. TANGENTIAL SHEAR PROFILE STACKED OVER FIVE 
SUBCLUMPS 

Since weak lensing signals of Coma cluster at the low 
redshift is weak and the number of background galax- 
ies within small radius is few, it is quite difficult to 
measure tangential shear profile for each subclump can- 
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Fig. 7. — The E- and B-mode maps reconstructed by background LSS model alone. Top-left : E-mode map in central region. Top-right : 
B-mode map in central region. Bottom-left : E-mode map in outskirts. Bottom-right: B-mode map in outskirts. The box sizes and contour 
scales are the same as those in Figure [2] and |3] Gray solid and black dashed lines show contours derived from LSS and observed shears, 
respectively. There are three LSS structures with a size of a few arcmins at a ~ 2ct level. In particular, the LSS convergences at a ~ Icr 
level are found in the clump candidates 5 and 8. 



didate. We therefore measure a profile of tangential 
shear components ensembling five subclump candidates. 
A mass measurement with stacked tangential distortion 
profile is complementary to the Cc statistics (^5]), because 
different shear catalogue is used in two measurements. 
In the statistics, source galaxies outside a given ra- 
dius is used, while, in the tangential shear measurement, 
source galaxies from inner to outer radius is indepen- 
dently available. We here exclude two dark halos as- 
sociated with cD galaxies in order to avoid a contami- 
nation of lensing distortion caused by the main cluster. 
The candidate 5 at which the projection effect is signif- 
icant is also ignored. The center for each subclump is 



chosen as the same position of the measurements 
(S|5|). The averaged tangential shear distortions of source 
galaxies, {g^){9n), is calculated in the circular annu- 
lus of the same radius, based on the same procedure as 
S|3l The typical projected distance between a center of 
a stacked tangential shear profile and the main cluster 
center, (Oos), is obtained with a weight function of lens- 
ing signals, (^off')' = E,(.9+,,)'^off,,/ E,(.9+j)' ^ 17^4, 
where (g-f,j) is the lensing signal (Sj3|) for each subclump 
and Ooffj is an angular separation between each sub- 
clump and the main center. 
We compute the LSS-corrected shear profile, 

{g+){9n) — {g'+^^'')(On), where the azimuthal aver- 
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TABLE 3 

Projected Mass for Subclump candidates 



ID 


S/N 


V 


^spur 


f'map 


Man 




(RA, Dec) 




Name 










% 




10i2/i-ij\^Q 


lOi2fc-iM0 


deg 






arcmin 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




(9) 


(10) 


1 


3.1 


3.6 


0.0 


0.0 


8.01 ± 1.05 


7.70 ± 1.05 


(195.042,27.996) 




NGC 4889 


1.2 


2 


3.9 


3.7 






4.28 ±0.83 


4.56 ±0.83 


(194.882,27.948) 




NGC 4874 


1.2 


3 


3.2 


2.9 


0.4 


1.3 


5.03 ± 0.87 


5.29 ± 0.87 


(194.980, 27.843) 




SA 1656-054 


1.2 


4 


3.6 


3.2 


0.0 


0.0 


3.49 ± 0.74 


3.73 ± 0.74 


(194.722, 27.829) 


SDSS 


J125848.72+274837.5 


0.7 


5 


3.4 


3.1 


12.2 


41.1 


5.95 ± 1.08 


2.74 ± 1.08 


(194.732, 27.722) 








6 


3.3 


3.1 


0.3 


1.0 


4.63 ± 0.84 


4.54 ±0.84 


(194.724, 27.612) 


SDSS 


J125858.10+273540.9 


1.4 


7 


3.1 


3.1 


0.1 


0.3 


6.07 ± 1.24 


7.32 ± 1.24 


(194.643, 27.625) 




NGC 4853 


1.7 


8 


3.3 


3.4 


1.8 


8.2 


5.90 ± 1.32 


4.61 ± 1.32 


(194.475, 27.044) 


SDSS 


J125756.65+270215.0 


0.8 



Note. — Col. (1): Name of subclump candidate. Col. (2): Signal- to- noise ratio in mass maps. Col. (3): v — i^/i^T-ms- The 
is obtained by 1000 Monte-Carlo simulations. Col. (4) Probability of spurious lensing peak considering contributions of main cluster, LSS 
effect and intrinsic noises. Col. (5): Bias by LSS effect and intrinsic noise for a preference to detect a clump in the convergence field of 



main cluster, defined as On 



3(LSS) , , p(noiso). 



7(PapurO- Col. (6): Projected mass for subclump without considering background LSS lensing 
effects. Col. (7): Projected mass for subclump with a correction of LSS lensing effects. Col. (8); Central positions for mass measurements. 
Col. (9); Name of galaxies in the mass clump candidate. Col. (10): Angular distance, d^ff, between luminous galaxy in clump candidate 
and centroid position for mass measurements. We note that they are consistent with galaxy positions within the uncertainty of centroid 
position for clump candidates (FWHM ~ 2f). 
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Fig. 8. — Aperture mass profile, Mq^, for each subclump candidate (1-8). Open and filled circles represent aperture masses without 
and with LSS correction, respectively. Aperture masses are saturated on the outer radii, which indicates that the mass density in the 
surroundings are quite low. It is consistent with results of stacked lensing analysis ([[B}. The aperture mass for candidate 5, corrected by 
LSS lensing effect, reduces to about half of that without a LSS correction. 
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age of the LSS distortion components, {g\_ )(^n), 
are calculated without statistical weight (ug.i = 1). 
Figure ^ shows the stacked tangential shear profiles as 
a function of transverse separation 6, with and without 
the LSS effect. We estimate the contribution of the main 
cluster mass on the stacked lensing signals, because 
the tangential shear provides full information on the 
lensing signals of gravitational potentials of both the 
main cluster and subclumps. It is necessary to calculate 
lensing distortions caused by the main cluster in order to 
measure typical mass of interior subhalos. We follow the 
convolution technique of Yang et al. (2006) to measure 
the azimuthally averaged convergence at the offset from 
the main cluster center in the lens plane. The values of 
g+ of the offset main halo is at the order of 0(— 10~^) in 
the range of the stacked lensing profile {9 < 6f), which 
indicates that the lensing signals from the main cluster 
at the positions of subclumps are negligible. 

Since the tidal field of main cluster disrupts dark mat- 
ter halos of interior substructures, the subhalo radius 
would be determined by the tidal radius rather than the 
virial radius rvh-- We therefore consider a truncated SIS 
model (TSIS) and a truncated NFW model (TNFW) for 
the tangential fitting, whose density profiles are trun- 
cated at the radius rt- 



Ptsis(?') ^ 



'2^ 

: for r > Tt 



M 



(TNFW) 



1 



(7) 



for r < Tt 



= for r > rt, (8) 



where cr^ j is a velocity dispersion for TSIS model, 

Mg™^^^ and Ct are a mass and a concentration for 
TNFW model, and r^,* — rt/ct is a scale radius de- 
termined by the concentration and the truncated ra- 
dius rt- The subclump mass for TSIS model is given 

Analytical expressions of the two-dimensional projec- 
tion S of the density field are obtained by an integration 



i(^)arccos(^) for 9 



< 



for 



> 



where 6 — r / Di is an angular size of the three dimen- 
sional radius, and 9E,t = 47r(o'i,.t/c)^£'is/£'s is the Ein- 
stein radius for TSIS model. 



V^t I 1 nrrrn-.Vi ( ^"^ 

(i-x2)(i+ct) + (i-x2)3/2arccohn i^^j^^^^j^ 
for a; < 1 



3(l+ct) l^-^ + 1+Ct ) 



for X ~\ 



(10) 



" (l-x2)(l+Ct) (2:2_l)3/2 



for \ < X < Ct 
for Ct < X 



The expression of /(x) is the same as Takada & Jain 
(2003) and Hamana, Takada & Yoshida (2004), although 
they adopt that the mass density is truncated at the 
virial radius. We here do not assume rt — r^„, because 
we aim to investigate disrupted interior substructures. 
The TSIS and TNFW profiles are specified in terms of 
two parameters of 9E,t a-nd 9t and three parameters of 

Mg™^^\ Ct and 9t, respectively. 

We fit the LSS- uncorrected and corrected distortion 
profiles with TSIS and TNFW models. The best-fit pa- 
rameters are summarized in Table SI Two models well 
describe the stacked tangential shear profile. The best- 
fit parameters with and without a LSS lensing correction 
does not change significantly. The truncated radii for two 
models are in good agreements with each other. Indeed, 
the break in the tangential shear profile is clearly found 
sX9 ^ If 75 (Figure ini). The values of 5+ steeply decrease 
(cx 0~^) over the truncated radius, which indicates that 
the halo mass density drops to zero at the truncated 
radius. It is consistent with the Cc mass measurement 
(^. This feature does not clearly appear in a stacked 
lens analysis for massive clusters of Mvh- > IO^/i-^Mq 
(Okabe et al. 2009). The TNFW and TSIS masses 
are in agreements with the the mean projected mass 
(M^d"^)34678 = (4.78 ± 1.07) X 10i2/i-iAfQ. We note 
that the two-dimensional mass for truncated mass mod- 
els (TNFW and TSIS) has the same analytical expres- 
sion as the three-dimensional mass (Afj'jf — A/g^j^^^ and 
Af™^^ — Af™^^), because there is no projection ef- 
fect due to the zero mass outside rj. 

We also investigate model fittings for spurious mass 
clumps appeared in mass maps simulated by rotated 
shear catalogue. First, we randomly rotate an angle 
in the ((71,52) plane for each background galaxies with 
it's \g\ fixed and then conduct the mass reconstruction 
for background catalogue by 500 times. Second, we de- 
tect peaks whose significance is over Sti in mass maps. 
Third, we measure stacked tangential profiles for 5 peaks 
which are bootstrap re-sampled from detected peaks by 
300 times. In the measurements, the peak whose dis- 
tance from cluster center is more than lOf is only used in 
order to avoid the main cluster lensing signal (Yang et al. 
2006). The averaged x^s are (x^) = 9.9 for TNFW and 
(x^) = 15.7 for TSIS models. Those x^s are worse than 
our results. This is why the form of stacked tangential 
profile for spurious clumps is different from Figure [9l 



7. DISCUSSION AND CONCLUSION 
, r(TNFW) ^, , , ^ ^ 

in\ -''^snh CI \ nin /n\ 7.1. Cluster Mass Comparison 

^TNFW^H o ^ 2 ^ ' x^^^s,t 9 

27rEcr?';.tm(ct) 
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Fig. 9. — Profiles of tangential shear component (top panel), g^, and the 45 degree rotated component (bottom panel), gx , obtained from 
an ensemble of five subclumps. Left : the stacked profile from original shear catalogue. Right : the profile with a LSS lensing correction. 
The 3+ values clearly decrease over 9 ~ lf75. The solid and dashed lines are subhalo models of TSIS and TNFW, respectively. 



TABLE 4 

The best-fit truncated SIS (TSIS) and truncated NEW 
(TNEW) MODEL parameters. 



TSIS 


LSS uncorrected 


LSS corrected 


(1) 


(2) 


(3) 


8t [arcmin] 
rt [/i~^kpc] 
(T„,t [km/s] 


o 47+1.21 
-0.56 

^- '"-0.15 
or; 91+11.32 
OO.Zl_2 98 

460.59l30:« 


1:4+1.32 

1 76+0-61 
1- '"-0.13 

or Qi +12.25 

oo.zi_2 58 

465.27^30.21 


xVd.o.f 


4.4/7 


5.3/7 


TNEW 


LSS uncorrected 


LSS corrected 


(4) 


(5) 


(6) 



[10-h-iMe] 




4 nn+1.45 


9t [arcmin] 


, 7c:+0.25 
^- ''-'-0.26 


, 7C.+0.21 
1. 'O-0.29 


rt [/i~lkpc] 


35.011^:11 


35.0llt.?^ 


Ct 


1.00_Q 79 


, rq+0.74 


xVd.o.f 


2.5/6 


2.4/6 



Note. — Col. (1): Name of best-fits parameters for TSIS 
model (eq. O : the subhalo mass ^'^^^^\ the truncated radius, 
rt, its angular radius dt and the velocity dispersion a-u^t- Col. (2-3) 
: The best-fit parameters for TSIS model without and with LSS 
correction. Col. (4) : Name of best-fits parameters for TNEW 

model (eq. (S)!: the subhalo mass ^guif^^'i the truncated radius, 
rt, its angular radius dt and the concentration parameter ct- Col. 
(5-6) : The best-fit values for TNEW model without and with 
LSS correction. All parameters are determined by fitting the mean 
distortion profile which is obtained by staking the distortion signals 
for 5 clump candidates (3,4,6,7 and 8). The x^/ti-o.f is the chi- 
square for best-fits and the degree-of-freedom. 



We compare our mass estimates with previous results 
of multi-wavelength data. The line-of-sight velocity dis- 
tribution of member galaxies with the Jeans equation 



which requires the assumption of the dynamical equilib- 
rium derived Mvir = 9.8 x lO^'^/i'^Af© and Cvir = 9.4 for 
NFW mass model (Lokas & JVIamon 2003). The caus- 
tic mass estimates to use a characteristic pattern in the 
redshift-space, formed by galaxies falling into cluster po- 
tential (Kaiser 1987) obtained M200 = 7 .d^^y^W^h'^ Mq 
and C200 = 10 (Rines et al. 2003). Our resuh of NFW 
mass Afvir and Af2oo is in good agreement with their re- 
sults, while our concentration parameter is lower. Our 
weak lensing analysis is to use one band imaging data. 
As demonstrated by Broadhurst et al. (2005), see also 
( Umetsu & Broadhurst 2008 ; Okabe et al. 2009), the 
dilution contamination of member galaxies on lensing sig- 
nals is problematic for an accurate measurement of the 
concentration parameter. It is therefore of importance 
to correct the dilution effect by secure selecting of back- 
ground galaxies in the color-magnitude plane (Broad- 
hurst et al. 2005; Umetsu & Broadhurst 2008; Umetsu et 
al. 2009). In addition, we require the data, which cov- 
ers the virial radius, in order to improve measurement 
accuracy of halo mass. The SDSS and CHFT weak lens- 
ing results (Kubo et al. 2007; Gavazzi et al. 2009) are 
M2nn = 18.8t^ixlO"MQ and M200 = 5.llt?x lO^Mfr 

+ 13.16 



and C200 = 3.84+|!)^;^g^ and C200 = 5.0_2.5 



-,+3.2 



0j 

They do 

not disagree with our results within large errors. The 
ASCA and ROSAT X-ray observations with assump- 
tions of hydrostatic equilibrium, isothermally and single 
(3 model shows M.soo = 11.99ti:29 ^ lO^^'M© (Reiprich 
& Bohringer 2002) and Afgoo = '^■^^tlH x lO^Af© 
(Chen et al. 2007), which are higher than our esti- 
mates. We cannot rule out a possibility that the low 
angular resolution of ASCA satellite leads to a bias on 
mass estimates. It is of critical importance to compare 
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X-ray and weak-lensing masses of Coma cluster which is 
the only cluster known to have turbulence in the ICM. 
The ASCA and XMM-Newton X-ray observations (e.g. 
Watanabe et al. 1999; Arnaud et al. 2001) have shown 
the complex temperature variations in the intra-cluster 
medium (ICM). Schuecker et al. (2004) has revealed a 
Kolmogorov/Oboukhov-type turbulence spectrum in the 
ICM as a consequence of the projected pressure distribu- 
tions. They constrained that the lower limit of turbulent 
pressure accounts for 10 percent of the total ICM pres- 
sure. Recent hydrodynamic N-body simulations pointed 
out that X-ray mass estimates with an assumption of 
hydrostatic equilibrium are biased low due to ICM tur- 
bulence, because the gas motion pressure of turbulence 
as well as bulk motions supports a part of the total pres- 
sure (e.g. Ervard et al. 1999; Nagai et al. 2007). Their 
results cast a doubt on accurate cluster mass measure- 
ment by X-ray analysis alone, which is a serious concern 
for cluster-based cosmological probes (e.g. Vikhlinin et 
al. 2009a,b; Okabe et al. 2010; Zhang et al. 2010). 
Therefore, a comparison of independent mass estimates 
is of great importance to understand, in a quantitative 
manner, how much the gas motion pressure affects the 
X-ray mass estimates. 

7.2. Comparison with X-ray image 

We compare an X-ray exposure-corrected image re- 
trieved from archival data of XMM-Newton with mass 
counters of the central region (Figure [T0|) . The XMM- 
Newton data in the outskirts region is lacked. The X-ray 
image shows some point sources associated with galaxies 
in clump candidates. The intra-cluster-medium (ICM) 
distributions are not correlated with mass clump can- 
didates. This is why the intra-cluster plasma can be 
escaped from the gravitational potential of subclumps. 
The sound velocity of the ICM is given by 



X— ray image (XMM-Newton) + mass contours 



1/2 



1457 



8.25kev) 



1/2 



km/s, (11) 



where the temperature ksT ~ 8.25 ± O.lOkeV (Arnaud, 
et al. 2001), the mean molecular weight fi = 0.62, and 
the proton mass rrip. The sound velocity is higher than 
the typical escape velocity from subclumps, as below 



Cs > Vc 



1/2 



991km/s. (12) 



Since the temperature of point sources is low ksT ^ 
IkeV and high metal abundance (Vikhlinin et al. 2001), 
it's sound velocity is lower than the escape velocity. 

7.3. Tidal Radius 

The stacked tangential shear profile for five subclumps, 
excluding subhalos associated with cD galaxies and can- 
didate 5, is well described by the TSIS and TNFW 
models. The fitting result gives the truncated radius 
Tt ~ 35/i~^ kpc, which coincides with galaxy-galaxy 
lensing results in clusters (e.g. Natarajan & Springel 
2004; Natarajan, De Lucia & Springel 2007). The trun- 
cated radius is much smaller than a truncated radius 
^ 200 kpc obtained by galaxy-galaxy lensing studies of 
fields (e.g. Hoekstra, Yee & Gladders 2004). It would 




RA [cleg] 

Fig. 10. — X-ray exposure-corrected image, no-subtracted back- 
grounds, overlaid with mass contours of the central 26.' X 26f cluster 
region. 



be due to the strong tidal field of the main cluster grav- 
itational potential. The tidal radius of a subhalo or- 
biting in a spherically, symmetric mass distribution of 
a cluster is obtained by the balance between the tidal 
force of the primary halo and the gravity of subclump, 
rtidai = (M3ub/(M(< rp){2-d\ogM/d\ogrp))^/^rp (e.g. 
Tormen et al. 1998), where rp is the pericenter radius 
which is the minimal radius from the cluster center dur- 
ing its orbiting history. Since we do not constrain its 
pericenter radius directly from a current position and 
not derive the three-dimensional radius, we instead as- 
sume the mean, projected offset radius (^off) of sub- 
clumps in stacked lensing analysis. Here, we assume the 
NFW and TNFW mass models for the main cluster and 
substructures, respectively. We obtain the tidal radius 
''tidal ~ 42 /i~^kpc, which coincides with the truncated 
radius rt ~ 35 /i~^kpc. 

7.4. Subhalo Mass Fraction 

We found four subhalos and two cD galaxy halos in the 
central data (r ;$ 30f) and one halo in the outskirts data 



(30.' , 



60.'). There is a difference of the halo number 



for the radius, which might support the results of numer- 
ical simulation (e.g De Lucia et al. 2004;Gao et al 2004) 
that the number of substructure increases as the radius 
decreases. Compensating the limitation of our data re- 
gion, we roughly estimate the number of subhalos within 
the radii rvir and r2oo ■ If we assume that there are four 
(r £ 30.') and one (30f £ r ^ 60'.) subhalos in the area 
corresponding to our data, the halo number is estimated 

to be 7r(r2i^-30.'^)/Aoutskirts + 4x (^30.''Mccntor) + 2cDs, 

where A is the area of our data. The Poisson noise for 
distributions is applied for the statistical errors. 

The halo number detected by weak lensing analysis 
is expected to be Nsuh{< r^ii) = 43 ± 30 and iVsub(< 
r2oo) = 27 ± 15. If the typical halo mass is the best-fit 

value (Aff™^^^) = 4.00x lO^^/i-iMo obtained from fit- 
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ting of the stacked tangential profile, the total substruc- 
ture mass within the virial radius account for '-^ 19 ± 13 
and 17 ± 9 percents of total cluster masses Mvir and 
-^200, respectively. Although the total mass fraction con- 
tained in subhalos does not agree each other among lit- 
erature (e.g. De Lucia et al. 2004; Natarajan, De Lucia 
& Springel 2007, Gao et al. 2004), most authors esti- 
mate 5 -20 percent. Our result is in rough agreement 
with the numerical simulations. A galaxy-galaxy lensing 
study (Natarajan, De Lucia & Springel 2007) indicates 
that 10 — 20% of the mass is contained in cluster sub- 
structures, which also roughly agrees with our result. 

Our weak lensing analysis on the nearby cluster would 
indicate the possibility that the mass function of clus- 
ter substructures is measurable without assumptions of 
mass-to-light ratio for member galaxies and dynami- 
cal state, while the galaxy-galaxy lensing studies (e.g. 
Natarajan & Springel 2004; Natarajan, De Lucia & 
Springel 2007) requires the assumption of the mass and 
light scaling law. We however have not yet obtained 
the mass spectrum as the galaxy-galaxy lensing studies 
(e.g. Natarajan & Springel 2004; Natarajan, De Lucia & 
Springel 2007). 

Alternative possible approach to investigate the mass 
function of cluster substructures is to measure higher or- 
der moments of the lensed images (HOLICs) and using 
the moments to estimate the flexion (e.g. Okura et al. 
2007; Okura & Futamase 2008; Okura et al. 2008). They 
have shown for the first time that flexion analysis can dis- 
cover substructures using the image of A1689 (Okura et 
al. 2007). 

As pointed out by Shaw et al. (2006), the median 
mass fraction is a increase function of the virial mass 
(/sub oc M°jj*'*^°*'®), because massive objects, which 
formed more recently than less massive objects, have 
less time to disrupt subhalos (Zenter et al. 2005). The 



statistical study of mass fraction of galaxy clusters is, 
therefore, one of good tests of ACDM and the hierarchi- 
cal clustering, as the concentration-mass relation of the 
NFW mass model (Okabe et al. 2009). Hence, further 
systematic study of mass fractions is required. 

The area of our current data is insufficient to derive the 
halo mass function as well as to measure cluster mass ac- 
curately. The next instrument of a prime focus camera 
of Subaru telescope, Hyper-Suprime-Cam, whose field- 
of-view is ~ 1.5deg^, will efficiently observe the nearby 
cluster and enables us to conduct weak lensing and flex- 
ion analyses. Our result using the Subaru/Suprime- 
Cam does guarantee that weak lensing analysis using 
Subaru/Suprime-Cam and Hyper-Suprime-Cam is capa- 
ble for almost X-ray clusters. 
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